Polydispersity stabilizes Biaxial Nematic Liquid Crystals 
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Inspired by the observations of a remarkably stable biaxial nematic phase [E.v.d. Pol et ai, Phys. 
Rev. Lett. 103, 258301 (2009)], we investigate the effect of size polydispersity on the phase behavior 
of a suspension of boardlike particles. By means of Onsager theory within the restricted orientation 
(Zwanzig) model we show that polydispersity induces a novel topology in the phase diagram, with 
two Landau tetracritical points in between which oblate uniaxial nematic order is favored over the 
expected prolate order. Additionally, this phenomenon causes the opening of a huge stable biaxiality 
regime in between uniaxial nematic and smectic states. 
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Since its first prediction back in the early 1970s [TH3], 
the biaxial nematic (Nb) phase has strongly attracted 
the interest of the liquid crystal (LC) community [4]. 
In contrast to the more common uniaxial nematic {Njj) 
phase, where cylindrical symmetry with respect to the 
nematic director determines optical uniaxiality, the Nb 
phase is characterized by an orientational order along 
three directors and consequently by the existence of two 
distinct optical axes. The prospect of inducing orien- 
tational ordering along three directions, while maintain- 
ing a nematic fluid- like mechanical behavior [5], renders 
biaxial nematics preeminent candidates for next gener- 
ation LC-based displays |Bj. Although experimental ev- 
idences of stable Nb phases were reported already 30 
years ago in lyotropic LCs 7 , in thermotropics this re- 
sult was achieved in systems of bent-core molecules only 
a few years ago [S]. Actually, when trying to experi- 
mentally reproduce an Nb phase, one often encounters 
practical problems related to its unambiguous identifica- 
tion [1] and to the presence of competing thermodynamic 
structures Stabilizing Nb states is therefore an 

open, challenging scientific problem with huge potential 
applications. Motivated by the exciting results of a recent 
experiment on a colloidal suspension |12] . we use here a 
mean-field theory to investigate the role played by size 
polydispersity on the stability of biaxial nematics in sys- 
tems of boardlike particles. We show that a difference in 
the particle volume of a binary mixture can favor oblate 
uniaxial orientational ordering over prolate, in sharp con- 
trast with the behavior of the pure systems. This phe- 
nomenon gives rise to a new phase diagram topology 
due to the appearance of two Landau tetracritical points, 
leading to a wider region of Nb stability. This feature is 
shown to hold also for a larger number of components, 
thus offering an explanation to the results of Ref. [H]. 
Finally, we argue that our findings could furnish a new 
way to look for biaxiality in thermotropic LCs. 

At low density in lyotropics, and at high temperature 
in thermotropics, the Nb phase appears as a crossover 
regime in between "rod-like" and "plate-like" behavior 



[5]. In fact, one can distinguish between the Njj phase 
developed by rods, in which particles align the longest 
axis along a common direction (uniaxial nematic prolate, 
iV+), and that developed by plates, in which particles 
align the shortest axis (uniaxial nematic oblate, iV_). A 
natural candidate system for developing an Nb phase 
is a binary mixture of rods and plates |13j : however, in 
most cases a demixing transition into two uniaxial ne- 
matic phases, i.e. N^ and N-, prevents its stabilization 
[lOl [TT] . Alternatively, a stable Nb state is expected 
in a system of particles with cuboid (i.e. rectangular 
parallelepiped) shape defined by the lengths of the prin- 
cipal axes L > W > T, as depicted in Fig. [TJa) [3]. 
In this case, it is convenient to introduce a shape pa- 
rameter V, defined hy y = ^ — By increasing the 
packing fraction and disregarding the possible stability 
of inhomogeneous phases, a system of cuboids undergoes 
an / N+ Nb sequence of phases \i v > 0, whereas 
an / — 7^ N_ — >• Nb sequence is found if < (/ stands 
for the isotropic phase) [I3|. A schematic representation 
of these nematic phases is given in Fig. [ijb)-(d). The 




FIG. 1. (a) Cuboidal particle with dimensions L x W x T. 
(b) Schematic representation of a system of freely rotating 
cuboids in the biaxial nematic phase Nb, (c) the uniaxial 
nematic prolate A'^^. and (d) the uniaxial nematic oblate A'^-. 
In this work the rotational degrees of freedom are discretized 
according to the Zwanzig model [18] . 
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case V = describes the optimal "brick" shape exactly in 
between "rod- like" and "plate- like" . In this case the Njj 
phase is suppressed and substituted by a second-order 
INb transition [T3]. 

The first experimental realization of the hard-cuboid 
model was found only recently in a colloidal suspension 
of boardlike mineral goethite particles [T^] . By produc- 
ing particles with shape parameter ~ 0.1 close to zero 
((L) X (W) X (T) = 254 X 83 X 28 nm^ and size polydis- 
persity of 20 — 25%), the authors were able to produce an 
Nb phase stable over a pressure range surprisingly much 
wider than to be expected from theory [51 115j and simu- 
lations '16] for particles whose shape parameter deviates 
even slightly from zero. Even more interestingly, the au- 
thors affirm that no Njj phase was observed, contrasting 
Ref. [13] . They suggest that a possible reason for this 
disagreement should be found in ingredients whose effects 
have never been studied so far because of their complex- 
ity, i.e. fractionation, sedimentation and polydispersity. 
These unexpected results motivate our interest in ana- 
lyzing the effect of the above mentioned ingredients, in 
particular polydispersity, on the stability of the Nb phase 
in a fluid of hard cuboids. 

We consider an M-component suspension of col- 
loidal cuboidal particles of species a = 1,...,M with 
dimensions x Wa x Tq, {La > Wa > Tq) in a 
volume V at temperature T. The total number den- 
sity of colloids is n = J^a ^a/V, the mole fraction of 
species a is = Na/{nV) and the packing fraction is 
r] = n^^XaLaWoTa- The theoretical framework used 
in this Letter consists of Onsager theory of LCs [T7j, 
which is a density functional theory truncated at second- 
virial order. In order to facilitate the calculations we 
follow Zwanzig by restricting the orientations of the par- 
ticles to the six in which their principal axes are aligned 
along a fixed Cartesian frame [18]. Although quantita- 
tive agreement with real systems is not expected because 
of the simplifications introduced in the model, the same 
model was shown to successfully predict non-trivial phe- 
nomena like demixing in rod-plate mixtures |10| , orienta- 
tional wetting due to confinement and capillary nemati- 
zation (19j . Moreover, we expect that transitions between 
different nematic phases and smectic phases are better 
described by this model than transitions from isotropic 
to nematics. In density functional theory the free energy 
of the system is expressed as a functional of the local 
density pf(r) of particles of species a — 1,...,M with 
orientation i — 1, 6 as [20] 



= / dvY^pt{v)\\n{pt{v)Kl)-l 
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where fc^ is the Boltzmann constant and the thermal 
volume of species a. At second- virial order the excess 
free energy T"^^ reads 



dr dr' 



E 



/^,"'(r-r')p?(r)p^ (r'), 



(2) 

where /^P'(r) = exp[-<?'(r)/(A:BT)] - 1 is the Mayer 
function, defined in terms of the pair-wise potential 
wfj" (r). By neglecting spatial modulations, i.e. by im- 
posing pf (r) = pf, the free energy Eq. ([T]) reduces to 
an Onsager-type functional whose minimization (under 
the constraints that J^iPf — "-^a for all a = 1,...,M) 
allows to identify the spatially homogeneous equilibrium 
phase (see Appendix A) . Since at sufficiently high density 
one expects spatially inhomogeneous phases to be ther- 
modynamically favored, we apply bifurcation theory [14] 
to determine the limit of stability of the homogeneous 
equilibrium phases with respect to smectic fluctuations. 
By considering spatial density modulations only along 
the z axis, i.e pf (r) = pfiz) in Eq. ([T]), the smectic bi- 
furcation density is the minimum density at which the 
Hessian second-derivative matrix of the free energy has 
an eigenvalue equal to zero (see Appendix B). 

Our analysis starts by considering the simplest case 
of polydispersity, i.e. a mixture of M = 2 components 
with mole fractions Xi and X2 — 1 — Xi, respectively. 
Among the different ways one can parameterize polydis- 
persity, our preliminary analysis suggests to consider vol- 
ume polydispersity (i.e. same particle shape but different 
volume). Therefore, we study the phase behavior of a 
binary mixture of hard cuboids whose dimensions are 



Li=L{l + s), Wi = W{l + s), Ti=T{l + s), 
L2=L{l-s), W2 = W{l-s), T2=r(l-s), 



(3) 



where the parameter s g [0, 1) describes the degree of 
bidispersity. Notice that Eq. (13]) implies the same as- 
pect ratios for both species Li/Ti ~ L2/T2 — L/T and 
Wi/Ti = W2/T2 = W/T (hence vi ^ V2 = v). Here we 
set L/T = 9.07 and W/T = 2.96 {u = 0.1) in order to 
reproduce the experimental system of Ref. [12] , thereby 
neglecting the small effect of the ionic double layer used 
by the authors to interpret the experimental data. 

Fig. [2] shows density-composition phase diagrams of 
binary mixtures (M = 2) of boardlike particles with the 
experimental shape parameter v = 0.1 for various bidis- 
persity parameters (a) s — 0.15, (b) 0.18, (c) 0.20 and 
(d) 0.30, featuring isotropic (J), uniaxial nematic (A^+ 
andN^), biaxial nematic [Nb) and smectic [Sm) phases. 
Due to the near-perfect "biaxial" shape of the particles, 
fractionation is extremely weak and invisible on the scale 
of Fig. [2] (see Appendix C). At the extreme mole frac- 
tions xi = and Xi = 1 (pure systems) all phase dia- 
grams feature the phase sequence / — ^ Sm that 
is well known and expected for board-shaped particles 
with ly > 0, with the Nb phase metastable with respect 
to the Sm phase [9] [14] (see also Appendix D). How- 
ever, for all s > there is an intermediate composition 
regime in which the Nb phase is found to be stable, the 
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FIG. 2. (Color online) Phase diagram of a binary mixture of 
hard cuboids in terms of packing fraction r\ vs. mole fraction 
of the larger component x\ showing isotropic (/), uniaxial 
(A'^+ and A'^-) and biaxial [Nb) nematic and smectic (Sm) 
phases. The size of the particles is defined by Eq. ([3| with 
LjT = 9.07, WjT = 2.96 (y = 0.1) and bidispersities (a) 
s = 0.15, (b) s = 0.18, (c) s = 0.20, (d) s = 0.30. The 
solid lines separate different homogeneous phases, the dashed 
lines indicate the limit of stability of the homogeneous phases 
with respect to smectic fluctuations, whereas the open circles 
represent the Landau tetracritical points. 



more so for increasing s. Whereas the opening-up of 
a stable Nb regime is only quantitative for s — 0.15, 
there is a qualitative change of the phase diagram topol- 
ogy beyond s = 0.18, where two Landau tetracritical 
points appear (open circles in Figs. [2]jb)-(d)). In be- 
tween these critical points a region of stable 7V_ phase, 
which is not expected for the rod-shaped particles {v > 0) 
of interest, opens up. Clearly, Figs. [2]jc) and (d) show 
that this unexpected N- regime enlarges with bidisper- 
sity, accompanying a further increased Nb stability. In 
other words, excluded-volume interactions in mixtures of 
board-shaped rods with the same shape and different vol- 
ume tend to favor Nb stability as a consequence of an 
unexpected N^ — N_ competition. At higher packing 
fractions the increased Nb stability with respect to the 
Sm phase is not a surprise, given that regular packing 
into layers is hindered by size differences between parti- 
cles [11]. 

It is interesting to understand how the remarkable fea- 
tures of the binary mixture described in Fig. [2] change 
with the shape of the particles. Here we are mainly in- 
terested in two properties of the phase diagram: (i) the 
minimum threshold bidispersity sthr at which the Landau 
tetracritical points appear and (ii) the tetracritical mole 
fractions xl in terms of the bidispersity s. We change the 
particle shape (i/ = L/W - W/T) by fixing in Eq. ^ 
one aspect ratio (W/T) and varying the remaining one 
(L/T). Fig. [3]; a) shows for W/T = 2.0, 2.96, 4.0 and 5.0 
a similar trend: the minimum threshold bidispersity Sthr 
increases the more the shape deviates from the optimal 
"brick" one. At the same time, the fact that at fixed ly 




3. (Color online) (a) Threshold bidispersity Sthr for the 
appearance of a tetracritical point as a function of the shape 
parameter i/ for different fixed values of W/T. (b) Critical 
mole fraction xl as a function of the bidispersity parameter 
s for a binary mixture of hard cuboids for different shape of 
the particles (cf. Eq. ((sl). 



the threshold bidispersity decreases with W/T, indicates 
that the appearance of the Landau tetracritical points is 
favored by an increasing aspect ratio of the particles, in 
qualitative agreement with Ref. |21| . Moreover, by fixing 
the aspect ratio W/T — 2.96, we can observe the tetra- 
critical mole fraction as a function of the bidispersity for 
different values of = 0.01, 0.1 and 0.25 in Fig. ^h). 
The closer the shape is to the optimal "brick" , the wider 
is the difference in value of the two tetracritical mole frac- 
tions xl and, consequently, the stability regime of the iV_ 
phase. Finally, we note that no critical composition is ob- 
served if the particles are closer to the "plate-like" shape, 
i.e. it i^i — 1/2 = ly < one finds the N- in between the / 
and Nb phases for every value of s and xi (not shown); 
the N^ phase does not occur in this case. 

In order to analyze proper polydispersity, and thus 
more realistically model the experimental system of Ref. 
|12j . we extend our phase-diagram calculations to a sys- 
tem with M — 21 components of cuboids. Inspired 
by our analysis of the binary mixture and by the ex- 
periments 132], we fix the aspect ratios of all species to 
La/To, = L/T = 9.07 and W^/T^ = W/T = 2.96, such 
that (i) all species have the same shape = v = 0-1 and 
(ii) the size of each species is completely determined by 
Ta ■ We consider Ta to be distributed according to a dis- 
cretized Gaussian function with average (T) = 28 nm and 
standard deviation <j{T), where a is the size polydisper- 
sity. In general the calculation of a (high-dimensional) 
phase diagram of a multi-component system is a daunt- 
ing task [22]. In this case, however, it is justified to ignore 
fractionation (see Appendix C), which reduces the prob- 
lem to minimizing the functional with respect to at 
fixed nxa- The resulting phase diagram in the density- 
polydispersity representation is shown in Fig. Qa), fea- 
turing again /, N^, N_, Nb and Sm equilibrium states 
and a tetracritical point at tr ~ 24%, which is surpris- 
ingly close to the size polydispersity in the experiments 
jl2| . The strikingly large stability regime of the Nb is 
caused by the reduced stability of Sm and N^ (cf. Fig. 
|4][b)), not unlike in the binary case. However, a direct 
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ing this interpretation scheme, it is intriguing to wonder 
whether it is possible to enhance the Nb stabihty by con- 
sidering two- or muhi-component mixtures of molecules 
with biaxial symmetry and different size. We hope our 
findings will stimulate further research in this direction. 

This work is financed by a NWO-VICI grant and is 
part of the research program of FOM, which is financially 
supported by NWO. 

Appendix A: Density functional theory 



FIG. 4. (Color online) (a) Phase diagram of M = 21 compo- 
nents of hard cuboids (packing fraction rj) with aspect ratios 
Lc/Tc = 9.07 and wS/Ta, = 2.96 (a = 1,...,M) and Gaus- 
sian distributed dimensions with polydispersity a (see text). 
The dashed line indicates the limit of stability of the homoge- 
neous phases with respect to smectic fluctuations. The dotted 
rectangles highlight (b) the absence of the Nb phase at poly- 
dispersity o < 4% due to the direct N+Sm phase transition 
and (c) the reentrant character of the N+Nb transition close 
to the tetracritical point (open circle). 



INb transition similar to that observed in Ref. [T? is not 
expected in this model due to the reentrant character of 
the N+Nb phase transition (cf. Fig. |4|^c)). 

In conclusion, by means of a mean-field theoretical ap- 
proach with discrete orientations we have shown that 
size polydispersity strongly affects the phase behavior 
of boardlike particles, driving the emergence of a novel 
topology of the phase diagram. This topology change 
is due to the appearance of Landau tetracritical points, 
which in turn is related to a competition between the 
prolate "rod-like" ordering typical of the pure compo- 
nents and the oblate "plate-like" purely induced by the 
mixing. In combination with the destabilization of the 
Sm phase, we can conclude that polydispersity dramati- 
cally increases the stability regime of the Nb phase. The 
usual stability limitations oi Nb phases, such as N+ — N^ 
demixing of rod-plate mixtures and ordering into smec- 
tics, are therefore overcome in the present system. Al- 
though this work focuses on a particular value of the 
particles dimensions, its predictions hold for a more gen- 
eral choice of the relevant parameters, as reported in Fig. 
[3| Moreover, we do not expect the homogeneous phase 
behavior to be crucially dependent on the form of the in- 
teraction (cuboidal), on the contrary it should be qualita- 
tively similar to other excluded- volume interactions with 
the same symmetry (e.g. spheroid, spheroplatelet). 

Finally, it is tempting to consider this work in the per- 
spective of stabilizing Nb thermotropic liquid crystals. 
In this case, the soft-core character of the inter-molecular 
interactions does not allow for a univocal definition of 
"shape", and van der Waals forces can significantly in- 
fiuence the phase diagram. Nonetheless, it is widely ac- 
cepted that hard-core models contain the essential phys- 
ical ingredients for a first-approximation description of 
the structure of a molecular or colloidal fluid E3] . Follow- 



In the present work the orientational degrees of free- 
dom of the particles are treated within the Zwanzig model 
|18) , hence a particular orientation can be identified with 
a number z = 1, 6 (cf. Tab. |l|. 
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TABLE L Enumeration of the orientational configurations of 
a hard cuboid within the Zwanzig model. Each configuration 
i is identified with the directions {x, y, z) along which the 
particle axes (L, W, T) are aligned. 

According to density functional theory it is possible to 
express the free energy of a system as a functional of the 
single-particle density pf(j) of particles with orientation 
i {i = 1,...,6) belonging to species a {a = 1,...,M) as 



kBT 
where for brevity 



(Al) 



M 



dr = dr. 
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The excess term J^'^'^[p] has in general a non-trivial de- 
pendence on pf (r). For short-range potentials it is al- 
ways possible to express J-'^^[p] as a virial series in the 
single-particle density. Therefore, by truncating the se- 
ries at second-virial order and thus disregarding higher- 
order contributions, one obtains 



dvdv' /° 



'(r-r'K(r)p^ (r'), 
(A2) 
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where the Mayer function Z^"" (r) is defined in terms of 
the pairwise interaction potential uf^ (r) as 



/r (r) = exp 



° (r) 



(A3) 



The single-particle density pf (r) is related to the number 
of particles through the normalization condition 



(A4) 



For hard cuboids the interaction potential, which ex- 
presses the impenetrability of the particles, is 



The single-particle density at equilibrium is the one 
which minimizes Eq. (A7) with the constraints of Eq. 



(A8) for all a = 1, A/, hence it is found by solving the 



Euler- Lagrange equation 



a:„nexp(^-iEa'..'i^.v'p"') 



(AlO) 



which is achieved by standard numerical (iterative) tech- 
niques. 

Appendix B: Nematic-Smectic bifurcation 



" (r) 



oo if\x\<iXr+X^) 
and \y\ < (Y^'^ +Y^?') 
and \z\ < {Zf + Z<^')- 

otherwise. 



(A5) 



According to the index notation defined in Tab. [T] the 
6-dimensional vectors X", Y" and Z" of species a are 
given in terms of the dimension of the particles by 



While studying the homogeneous equilibrium phases of 
the system, we are also interested in estimating their up- 
per bound in the phase diagram, where spatially inhomo- 
geneous phases tend to be thermodynamically favored. 
Bifurcation theory US] provides a way to investigate 
the limit of stability of a particular phase. 

The condition of thermodynamic stability of a phase 
described by the single-particle density pf(r) requires 
that the system corresponds to a minimum of the free 
energy J^, i.e. a stationary point that satisfies 



X" = |(_Lq, Wq,, Tq, Lq, Wq,, Tq), 

Y" = ^(VKq, Tq, Lq,, Tq, £aj W^a)i 
Z" — \{Ton La, Wa, Wa, Ta, La). 



(A6) 



The main goal of this work is to study the stability 
of spatially homogeneous phases (i.e. isotropic and ne- 
matic). In order to simplify the problem we therefore 
neglect spatial modulations in the single-particle density 
by imposing pf{r) = pf. Consequently, Eq. (All be- 
comes 



VkeT ^ 

OL.% 



- 1 



E 

a, a' ,2,2' 



^tt' Pt Pi' > 



(A7) 

which is the restricted orientation version of th e Onsager 
free energy [T7] . The matrix Ef^ in Eq. ( A7 1 is the ex- 



cluded volume between two particles belonging to species 
a and a' with orientations i and i' interacting through 
the potential Eq. ( A5 ) 



drdv' 



Spfir)Spf,'{r') 



Spf{r)dpf, (r') > 0, 



(Bl) 

for any arbitrary perturbation Spf(r). By inserting the 
functional expression Eq. (lAll) into Eq. (Bl), one finds 



that the reference phase (described by pf{r)) ceases to 
be stable at the smallest density n = N/V at which a 
perturbation Sp°'{r) exists such that 



Spfir)=pnr) I dr' V/^,"'(r-r')5p?'(r'). (B2) 



Here we are interested in calculating the limit of stabil- 
ity of the (uniaxial or biaxial) nematic phase with respect 
to smectic fluctuations. With this in mind, in Eq. (B2) 



we neglect spatial modulations in the reference phase, i.e. 
pf{r) — pf, and a positional dependence of the fluctu- 
ations only along the z dire ction , i.e. Spf{r) — 6pf{z). 
After some rearranging Eq. (B2| becomes 



8(Xf + X? ) + Y^? ){Z^ + Z?). (A8) 



In the homogeneous case the normalization condition Eq. 
(A4) becomes 



E 



dz'QS? {z-z')a? (z'). 



where erf (z) — 5pf{z)/y/~f^ and 



(B3) 



(A9) 



Q?." {z)^Jpfpt: / dxdyf^?{v) 



(B4) 
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a symmetric (Hermitean) kernel. By inserting the ex- 
plicit form of the inter-particle potential (cf. Eq. (A3) 
and (|A5[)) into Eq. (|B4[), we obtain 



Q?/W=| if|z|<(Zf + 

1^0 otherwise. 

(B5) 

Eq. ( B3 1 can be more conveniently solved in Fourier 



space, where it reads 



(B6) 



with 



(?) - -ypfp^'^.v 'jo(g(^r + ^?')), (B7) 



and jo(a;) = sin(x)/a;. 

In conclusion, the limit of stability of the nematic 
phase with respect to smectic fluctuations can be numer- 
ically found as the minimum packing fraction rj* at which 
there exists a wave vector q* such that the 6M x 6M ma- 
trix with entries Q^? {q*) has a unit eigenvalue. The pe- 
riodicity of the corresponding bifurcating smectic phase 
is given by A* = 27r/g* . 
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FIG. 5. (a) Phase diagram of a binary mixture of hard 
cuboids in terms of the reduced osmotic pressure P* = 
PLWT /(kBT) vs. mole fraction of the first species x^. (b)- 
(d) / (red solid line), A'^+ (blue dotted line) and (green 
dashed line) branches of the Gibbs free energy per particle 
g = G/N at (b) P* = 0.6, (c) P* = 1.1 and (d) P* = 1.5. 
A straight line with slope a = dg / dxi_\xi=xi = dg / dxi_\x-^=xi^ 
(with xi and xn the compositions of the coexisting isotropic 
and nematic phases) was subtracted in each case to enhance 
the visualization of the common tangent construction. 



Appendix C: Nearly second-order character of the 
IN± transition 

When dealing with mixtures, the phase diagram is con- 
veniently expressed in terms of the osmotic pressure P 
vs. the mole fraction of M — 1 components. In this 
way it is possible to visualize the coexistence of phases 
characterized by a different composition with respect to 
the parent distribution. This phenomenon, called demix- 
ing or fractionation, is a consequence of the first-order 
character of the transition. 

Here we analyze demixing in a binary mixture of 
cuboids parameterized as in Eq. (3) with L/T = 9.07, 
W/T = 2.96 and s = 0.2. In Fig. [sja) we report the 
phase diagram for such a system as a function of the mole 
fraction xi of the larger species. The expected first-order 
character of the INjj transitions is not detectable at this 
scale (see below), whereas the NjjNb transitions appear 
to be second order. At three different values of the re- 
duced pressure P* = PLWT/ (kBT) we calculated the 
isotropic and uniaxial nematic branches of the Gibbs free 
energy per particle g{P,xi) = G{P, Ni, N2) / {kBT{Ni + 
7V2))- The coexistence between the two phases is given 
by a common tangent construction, which allows to eval- 
uate the difference in composition Axi of the coexisting 
phases. The results are reported in Fig. [5]jb)-(d) for 



P* = 0.6, 1.1 and 1.5, respectively. In the three cases, 
two of which describe a INj^ and one a /iV__ transition, 
Azi K, 10~^ and can therefore be neglected. The situa- 
tion does not change when one considers different values 
of the bidispersity parameter s. 

Although Landau-de Gennes theory predicts the INu 
transition to be first order |26| . we have just shown that 
its discontinuous character can be safely neglected for the 
binary mixture of boardlike particles we consider in this 
work. In our opinion, this fact is tightly related to the 
shape of the particles close to the v = Q value. In fact, 
when considering a monodisperse system, the closer v is 
to zero the weaker is the first-order character of the INjj 
transition (see also Sec. [d|. This fact allows us to assume 
that for an arbitrary number of components of volume- 
polydisperse cuboids close to v — Q the INjj transition 
can be approximated as continuous. As a consequence, 
we can neglect demixing in the phase behavior analysis 
reported in Fig. 4, thus reducing enormously the com- 
plexity of the problem. 
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FIG. 6. Phase diagram of a monodisperse system of hard 
cuboids as a function of the shape parameter v = LjW — 
W/T (with L/T = 9.07 fixed and W/T variable). The solid 
lines indicate phase boundaries as calculated by minimizing 
the Onsager-Zwanzig functional, the dashed line indicates the 
limit of stability of the nematic with respect to the smectic 
phase and the open circle the Landau tetracritical point. The 
inset highlights the first order character of the INu transition 
and how this tends to become continuous by approaching v = 
0. 



FIG. 7. Orientation distribution function of a monodisperse 
system of hard cuboids as a function of the packing fraction 77 
obtained by minimization of the Onsager-Zwanzig functional 
Eq. ( [A7| for M = 1. The cuboids have dimensions L/T = 
9.07 and (a) W/T = 3.04, (b) W/T = 3.01, (c) W/T = 2.99, 
(d) W/T — 2.96. The different lines indicate the probability 
of a particular orientation i — 1, 6 (cf. Tab. [l]|. The dashed 
vertical line shows the limit of stability of the nematic phases 
with respect to the smectic. 



Appendix D: Monodisperse system of hard cuboids 

The main goal of the present work is to investigate how 
polydispersity affects the phase behavior, and in partic- 
ular the stability, of the Nb phase in a system of hard 
cuboids. For this reason, it is instructive to study what 
the theoretical framework described in Sec. |X] predicts 
in the monodisperse case M = 1. In particular, we will 
focus here on the role of the particles dimensions on the 
phase behavior of the system. 

In Fig. [6] we report the phase diagram of a monodis- 
perse system of hard cuboids as a function of the as- 
pect ratio W/T at fixed L/T = 9.07. Consequently, 
by varying W/T one varies the shape parameter ly — 
L/W — W/T, in such a way that by crossing the point 
1/ = one expects a transition from plate- to rod-like 
behavior. This is precisely what Fig. |6] shows, where 
the phase separation lines are calculated by minimizing 
the Onsager-Zwanzig functional Eq. ( |A7[ ) with the con- 
straint of Eq. ( A9 1 for each value of the packing fraction 
77. Moreover, bifurcation theory (cf. Sec. |b]) provides a 
way to estimate the upper limit of stability of homoge- 
neous phases with respect to the smectic (dashed line in 
Fig. |6]). Fig. [6] shows that to observe a stable Nb phase, 
the shape of the particles should be designed with ex- 
tremely high precision in a small i^-regime about v — 0. 
In fact, for L = 9.07 T the Nb phase disappears unless 
2.96 T <W < 3.08 T. This is due both to the tight cusp- 
like shape of the NjjNb transition line and to the pre- 
empting character of inhomogeneous phases. Analogous 
results can be obtained by varying the shape parame- 



ter through L/T, while keeping W/T fixed (not shown). 
Finally, in the inset of Fig. [6] (note the different scale) 
we show the first order character of the INij transition, 
which tends to become second-order by approaching the 
critical point at = 0. 

For the sake of completeness, in Fig. [7] we report the 
orientation distribution function pi, which is the proba- 
bility of a given orientation i = 1, ...,6 as a function of 
the packing fraction rj for different values of the shape 
parameter v. In the monodisperse case this function co- 
incides with the single-particle density divided by the 
number density: pi = pi/n. The values of the orien- 
tation distribution function characterize the symmetry 
of the corresponding phase. In fact, at a given packing 
fraction 77 in Fig. [7] one can have one of the following 
possibilities: 

• the probabilities pi are all the same, i.e. pi = 1/6 
(isotropic / phase); 

• the probabilities Pi are coupled two-by-two, demon- 
strating the presence of a symmetry axis (uniaxial 
nematic Nu phase); 

• the probabilities Pi are different between each oth- 
ers (biaxial nematic Nb phase). 

Moreover, in the uniaxial nematic case one can further 
distinguish two situations: 

A the two more probable orientations have the short- 
est axis aligned along the same direction (uniaxial 
nematic oblate N- phase); 
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▲ the two more probable orientations have the longest 
axis aligned along the same direction (uniaxial ne- 
matic prolate 7V_|_ phase). 

This classification is easily generalized to the multi- 
component case. With this in mind, one can observe the 
difference in the orientation distribution function when 
V = -0.06 < {W/T = 3.04, Fig. ^a)), 1^ = 



(W/T = 3.01, Fig. ^h)) and v = 0.04 > {W/T = 2.99, 
Fig. [tJ^c)). The vertical dashed line indicates the limit of 
stability with respect to smectic fluctuations as given by 
bifurcation theory. Finally, Fig. shows the predicted 
orientation distribution function when the experimental 
value W/T — 2.96 is considered [12], and highlights how 
according to the model the Nb phase is expected to be 
preempted by inhomogeneous phases. 
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